Comparison between the ages for the genders (after controlling for type)

Questions

  • What are the differences between the ages for the different genders?
  • Do we observe the same changes as globally?

Age effect - General Questions

  • What are the differences between the ages?
  • Which genes and pathways are differentially expressed between 8w and 52w, between 52w and 104w, between 8w and 104w? Are they the same? Is there a gradient?
  • Are they different for the two genders?
  • Are they different for the two types?

Loads

Libraries and functions

Allowing multi-threading with up to 4 threads.
[1] "preparing gene to GO mapping data..."
[1] "preparing IC data..."
[1] "preparing gene to GO mapping data..."
[1] "preparing IC data..."
[1] "preparing gene to GO mapping data..."
[1] "preparing IC data..."

Data

  1. 14447
  2. 88
  1. 0.3609065973726
  2. 0
  1. 4.84577258072278
  2. 0.21654395842356

Differentially expressed genes

Warning message in pcls(G):
“initial point very close to some inequality constraints”Warning message in pcls(G):
“initial point very close to some inequality constraints”Warning message in pcls(G):
“initial point very close to some inequality constraints”Warning message in pcls(G):
“initial point very close to some inequality constraints”Warning message in pcls(G):
“initial point very close to some inequality constraints”Warning message in pcls(G):
“initial point very close to some inequality constraints”Warning message in stack.default(getgo(rownames(l$sign_fc_deg), "mm10", "geneSymbol")):
“non-vector elements will be ignored”Warning message in stack.default(getgo(rownames(as.data.frame(l$deg)), "mm10", "geneSymbol", :
“non-vector elements will be ignored”

Stats

All DEG (Wald padj < 0.05)All over-expressed genes (Wald padj < 0.05 & FC > 0)All under-expressed genes (Wald padj < 0.05 & FC < 0)DEG (Wald padj < 0.05 & abs(FC) >= 1.5)Over-expressed genes (Wald padj < 0.05 & FC >= 1.5)Under-expressed genes (Wald padj < 0.05 & FC <= -1.5)
52w VS 8w (F)21151190 925 881625 256
52w VS 8w (M) 393 252 141 248182 66
104w VS 52w (F) 239 146 93 159115 44
104w VS 52w (M)3848201718311420394 1026
104w VS 8w (F)17151020 695 915696 219
104w VS 8w (M)4189220619831884762 1122

All DEG (Wald padj < 0.05)

Some explanation (specially for the gender difference in microglia aging)

DEG (Wald padj < 0.05 & abs(FC) > 1.5)

DEG (Wald padj < 0.05 & abs(FC) > 1.5)

Log2FC

Z-score

Co-expression (WGCNA)

DEG into gene co-expression network

  • White: up-regulated
  • Black: down-regulated
Comp Male Female
52w vs 8w
104w vs 52w
104w vs 8w

Z-score in modules

GO analysis

Biological process

Dot-plot with the most over-represented BP GO (20 most significant p-values for the different comparison)

Using term, id as id variables
Using term, id as id variables

Network based on description similarity

Comp Male Female
52w vs 8w
104w vs 52w
104w vs 8w

52w VS 8w (F)

GO Tree at "../results/dge/age-effect/age_gender/go/52w_VS_8w_F.png"

52w VS 8w (M)

GO Tree at "../results/dge/age-effect/age_gender/go/52w_VS_8w_M.png"

104w VS 52w (F)

GO Tree at "../results/dge/age-effect/age_gender/go/104w_VS_52w_F.png"

104w VS 52w (M)

GO Tree at "../results/dge/age-effect/age_gender/go/104w_VS_52w_M.png"

Cellular components

Dot-plot with the most over-represented CC GO (20 most significant p-values for the different comparison)

Using term, id as id variables
Using term, id as id variables
Comp Female Male
52w vs 8w
104w vs 52w
104w vs 8w

Molecular functions

Dot-plot with the most over-represented MF GO (20 most significant p-values for the different comparison)

Using term, id as id variables
Using term, id as id variables
Comp Female Male
52w vs 8w
104w vs 52w
104w vs 8w

KEGG pathways

Error in `$<-.data.frame`(`*tmp*`, labels, value = c("", "", "", "", "", : replacement has 12 rows, data has 28
Traceback:

1. plot_kegg_pathways(age_gender_deg$over_represented_KEGG[, "category"], 
 .     age_gender_deg$fc_deg, "../results/dge/age-effect/age_gender/kegg/over_repr_kegg/")
2. suppressMessages(pathview(gene.data = fc_deg, pathway.id = cat, 
 .     species = "Mus musculus", gene.idtype = "Symbol"))
3. withCallingHandlers(expr, message = function(c) invokeRestart("muffleMessage"))
4. pathview(gene.data = fc_deg, pathway.id = cat, species = "Mus musculus", 
 .     gene.idtype = "Symbol")
5. `$<-`(`*tmp*`, labels, value = c("", "", "", "", "", "", "", 
 . "", "", "", "", ""))
6. `$<-.data.frame`(`*tmp*`, labels, value = c("", "", "", "", "", 
 . "", "", "", "", "", "", ""))
7. stop(sprintf(ngettext(N, "replacement has %d row, data has %d", 
 .     "replacement has %d rows, data has %d"), N, nrows), domain = NA)

Gender differences in aging

Question: Is there differences in aging between gender? Is there really a delay for some genes in male?

52w vs 8w for F 52w vs 8w for M 104w vs 52w for F 104w vs 52w for M Gene number
Set 1 != == == != 729
Set 2 == != != == 9
Gene number
Set 1676
Set 2 3

Genes with differential expression delayed in male

Genes (set 1):

  • Differentially expressed between 52w and 8w in F
  • Not differentially expressed between 104w and 52w in F
  • Not differentially expressed between 52w and 8w in M
  • Differentially expressed between 104w and 52w in M
  • No criteria for 104w vs 8w in both F and M

52w VS 8w (F) > 0 and 104w VS 52w (M) > 0

[1] 232

	Pearson's product-moment correlation

data:  mat[, 1] and mat[, 2]
t = 32.072, df = 230, p-value < 2.2e-16
alternative hypothesis: true correlation is not equal to 0
95 percent confidence interval:
 0.8773837 0.9251049
sample estimates:
      cor 
0.9040226

52w VS 8w (F) > 0 and 104w VS 52w (M) < 0

[1] 94

	Pearson's product-moment correlation

data:  mat[, 1] and mat[, 2]
t = -17.071, df = 92, p-value < 2.2e-16
alternative hypothesis: true correlation is not equal to 0
95 percent confidence interval:
 -0.9131327 -0.8127703
sample estimates:
       cor 
-0.8718158

52w VS 8w (F) < 0 and 104w VS 52w (M) > 0

[1] 181

	Pearson's product-moment correlation

data:  mat[, 1] and mat[, 2]
t = -9.1249, df = 179, p-value < 2.2e-16
alternative hypothesis: true correlation is not equal to 0
95 percent confidence interval:
 -0.6554453 -0.4549900
sample estimates:
       cor 
-0.5634548

52w VS 8w (F) < 0 and 104w VS 52w (M) < 0

[1] 169

	Pearson's product-moment correlation

data:  mat[, 1] and mat[, 2]
t = 18.845, df = 167, p-value < 2.2e-16
alternative hypothesis: true correlation is not equal to 0
95 percent confidence interval:
 0.7695674 0.8676545
sample estimates:
      cor 
0.8247168

Summary

Comp 52w VS 8w (F) < 0 52w VS 8w (F) > 0
104w VS 52w (M) > 0
104w VS 52w (M) < 0
149
179
124
56

Genes with differential expression delayed in female

Genes (set 2):

  • Not differentially expressed between 52w and 8w in F
  • Differentially expressed between 104w and 52w in F
  • Differentially expressed between 52w and 8w in M
  • Not differentially expressed between 104w and 52w in M
  • No criteria for 104w vs 8w in both F and M
1